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(1)  good  agreement  between  theory  and  experiment  with  regard  to  the  thrust  deduction 
of  a centerline  rudder.  (2)  at  equal  thrust  the  forward  and  aft  propellers  produced  73 
percent  and  27  percent  of  the  total  thrust  deduction,  respectively,  and  (3)  the  total  thrust 
deduction  is  reduced  by  unbalancing  the  propelling  thrust  with  smaller  thrust  carried  on 
the  forward  propeller. 
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NOTATION 


Propeller  disk  area 

D 

Drag  augmentation  coefficient  

2?rp  V2  R2 
T 

I Imist  loading  coefficient  

1/2  P(A0),  V2 

Thrust  loading  coefficient  

1/2  P(A0),  V2 

Augmented  drag  or  propeller  diameter 
Maximum  diameter  of  a body  of  revolution 
Axial  distance  between  CR  propeller  planes 
Froude  number  V/^/gL 
Force  (source) 

Circumferential  average  factor  for  interference  velocities 

Total  nondimcnsional  circulation  at  each  radius  of  one 
propeller,  ZF/ffDV 

Lerbs’  distance  factor 

Body  length 

Reference  length  (distance) 

Output  of  a point  source 

Output  or  strength  density  of  a surface  source 

Propeller  radius 

Directed  distance  in  genera! 

Surface  area 
Total  thrust  CR  set 
Propeller  thrust;  ratio  T,/T  and  T,/T 
Thrust  deduction  coefficient 
Axial  induced  velocity 
Tangential  induced  velocity 
Free-stream  velocity 
Propeller  speed  of  advanced  - w)V 
Local  axial  velocity  component 


'V 

Volume  mean  veloeity 

V 

C oaxial  veloeity  induced  by  point  sources 

w 

Wake  fraction 

X 

Propeller  radius  fraction 

5 

Average  contraction  ratio  at  aft, propeller 

6 

Influence  coefficient 

r 

Average  circulation  factor  at  af  propeller 

p 

Mass  density 

SUBSCRIPTS 

0 

Particular 

1 

Forward  prope'ler 

2 

Aft  propeller 

a 

Axial 

r 

Frictional 

ll 

Hub 

i 

Interference  induced 

P 

Potential 

s 

Self  induced 

l* 

Tangential 

ABSTRACT 


A theoretical  method  is  presented  for  calculating  the  steady  propulsive 
interaction  (thrust  deduction)  force  in  contrarotating  propellers.  Contrarotating 
propellers  operating  at  off-design  loading  and  spacing  as  well  as  the  contribution 
of  a rudder  were  investigated.  The  importance  of  the  separate  thrust  deduction 
of  the  forward  and  aft  propellers  in  analyzing  the  behavior  of  a CR  propeller  set 
was  shown.  Numerical  results  are  given  for  a MAR  AD  high-speed  containership. 
Some  principal  findings  for  the  subject  ship  arc:  (1 ) good  agreement  between 
theory  and  experiment  with  regard  to  the  thrust  deduction  of  a centerline  rudder, 
(2)at  equal  thrust  the  forward  and  aft  propellers  produced  73  percent  and 
27  percent  of  the  total  thrust  deduction,  respectively,  and  (3)  the  total  thrust 
deduction  is  reduced  by  unbalancing  the  propelling  thrust  with  smaller  thrust 
carried  on  the  forward  propeller. 
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Project  F 43  432,  Task  14438,  Work  Unit  1-1544-256. 


INTRODUCTION 

To  the  author’s  knowledge  the  important  and  intriguing  problem  of  the  drag  augmenta- 
tion (thrust  deduction)  of  a ship  with  contrarotating  propellers  has  not  been  solved  previously 
in  sufficient  detail  for  practical  design  purposes.1  The  problem  is  inherently  more  complex 
than  single-screw  propulsion  which  has  received  considerable  attention,2  10  mainly  because 


'llickling,  R„  “Propellers  in  lhe  Wake  of  an  Axisymmclric  Body,"  Institution  of  Naval  Architects,  Quarterly  Transactions, 
Vol.  99,  No.  4 (Oct  1957).  A complete  listing  of  references  is  given  on  page  29. 

2Tsakonas,  S.  and  W.R.  Jacobs,  “Analytical  Study  of  the  Thrust  Deduction  of  a Single-Screw  Thin  Ship,”  Stevens  Insti- 
tute of  Technology,  Davidson  Laboratory  Report  816  (Mar  1962). 

^Beveridge.  J.L..  “Thrust  Deduction  Due  to  a Propeller  Behind  a Hydrofoil,"  David  Taylor  Model  Basin  Report  1603 
(Oct  1962). 

4 

Beveridge,  J.L.,  “Effect  of  Axial  Position  of  Propeller  on  the  Propulsion  Characteristics  of  a Submerged  Body  of  Revo- 
lution," David  Taylor  Model  Basin  Report  1456  (Mar  1963). 

^Nowacki,  H..  "Potential  Wake  and  Thrust  Deduction  Calculations  for  Ship-Like  Bodies."  Transactions  STG  (1963)  (in' 
German). 

^Wald,  0.,  "Performance  of  a Proj>ellcr  in  a Wake  and  the  Interaction  of  Propeller  and  Hull,"  Journal  of  Ship  Research, 
Voi.  9.  No.  I (Jun  1965). 

Amtsberg,  II.,  "Investigations  on  the  Interaction  between  Hull  and  P--jpcller  of  Bodies  of  Revolution,"  David  Taylor 
Model  Basin  Translation  309  (Dec  1965). 

r 

Dregcr,  W.,  "A  Method  of  Calculation  of  Potential  Thrust  IX’duction,”  David  Taylor  Model  Basin  Translation  328 
(Mar  1966). 

a 

Pohl,  K.H.,  "The  Interaction  between  Hull  and  Propeller,"  David  Taylor  Model  Basin  Translation  334  (Tcb  1967). 

*°Bcvcridge  J.L.,  "Analytical  Prediction  of  Thrust  Deduction  for  Submersiblcs  and  Surface  Ships,"  Journal  of  Ship 
Research,  Vol.  13,  No.  4 (Dec  1969). 
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of  the  .uktition.il  singularities.  imolvcd  .uul  the  attemhint  multiple  integrations  required  for 
the  potential  par t of  the  thrust  deduction  force.  The  frictional  component  of  the  thrust 
deduction  is  known  to  be  very  small  and  can  be  approximated.11  13  This  elucidation  of 
the  problem,  i.e..  as  a l.agaltv  type  ol  force,  has  been  adequately  investigated  and  applied 
successfully  .3* 10  It  has  its  basis  in  a well  known  principle  of  reciprocity  (Newton)  for  finding 
the  forces  between  bodies  in  a How  field. 

In  apportioning  the  thrust  load  between  a pair  of  contrarotating  propellers,  current 
design  methods  do  not  consider  the  actual  division  of  the  thrust  deduction  between  the  forward 
and  after  propellers.  Thus,  the  design  condition  of  either  thrust  or  torque  balance  between 
propellers  at  equal  rates  of  rotation  (tpm)  may  not  produce  the  minimum  shaft  horsepower. 

The  investigation  rc  ported  herein  develops  a theoretical  method  for  calculating  the  thrust 
deduction  due  to  contrarotating  stern  propellers.  The  scope  of  the  work  includes  the 
development  of  analytical  methods  for  determining  the  effect  of  rudder,  variations  in  the 
ratio  ol  thrust  between  forward  and  after  propellers,  and  the  spacing  between  propellets. 
Computational  results  are  presented  for  a MAKAI)  high-speed  eontainership.14 

The  plan  and  prmc  pal  points  of  this  report  include:  (1 ) the  derivation  of  analytical 
expressions  for  the  »hrust  deduction  for  contrarotating  propellers,  and  (2)  the  introduction 
of  an  influence  coefficient  to  account  for  the  effect  of  wave-making  and  certain  geometric 
particulars  on  the  potential  wake  of  a high-speed  merchant-ship  type  of  hull. 

STATEMENT  OF  PROBLEM 

In  potential  How.  mathematical  singularities  which  represent  solid  bodies  or  boundaries 
are  widely  used.  The  Lagallv  steady -motion  equation15,16  relates  the  force  on  a source  (or 
sink)  of  given  strength  in  an  arbitrary  How  to  the  velocity  which  the  How  would  possess 
at  the  location  of  the  source  it  the  source  did  not  exist.  The  thrust-deduction  problem  is 
analogous  to  the  familiar  problem  of  finding  the  forces  between  bodies  in  a tlow  field.  A 
torce  arises  from  the  mutual  influence  of  the  hull  on  the  contrarotating  propellers  and  these 


1 1 Diskmann.  II. I . "I lie  Interaction  l«ctuecn  Propeller  .iml  Slop  with  S|>cu.il  Consideration  to  the  Influence  of  Waves." 
J.ilirliiK.li  tier  Silntfl.aiiteilmisiiicn  (•evvelvclull.  40  (1939). 

1 “Van  l-jinmcfcn.  VV.I*  A..  “ An.ilv.iv  of  Propulsion  Components  in  Kclulion  in  Si  ale  I 'fleet  by  Model  rests.”  David 
I a) lor  Model  llasm  translation  6K  (Sep  1950). 

1 3|ltveridce  I I ..  "Pressure  Iristnliiilioii  on  lowed  and  Propelled  Slreainline  llodics  ol  Revolution  al  Deep  Submergence." 
David  lav  lor  Model  llasin  Report  |fif>‘  (Jim  |96f>). 

14Slinm-  Icjscn,  J . ”A(  ompanson  of  Conir.iioialim'  Propellers  will)  oilier  Propulsion  Systems.”  Marine  Tcihnolojiy. 

V..I  9.  V.  1 ( I in  1972) 

^Durand,  VV  I . “ \ciodsn.imii  lliemv."  Div  (,  Vol.  I,  p.  2W).  Dover  Publications  (1963). 

16|fel/.  A . Dm  Mrlhoil  of  Sin-'iilarities  lor  ills'  IX'lenmnall.m  of  I ones  and  Moments  Ac  liny:  on  a llody  in  Potential 
I low."  David  I as  lor  Model  Hum  I r.msl.ition  ’41.  Revised  I dilion  (lun  1951). 
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propellers  on  the  hull.  Simply  stated,  a hull  a fleets  contrarotating  propeller  performance 
through  its  wake,  and  a hull  (from  the  point  of  view  of  resistance)  is  affected  by  the  self- 
induced  velocity  Held  of  the  contrarotating  propellers. 

Hie  Lagally  theorem  provides  a means  of  circumventing  the  detail  of  integrating  pressures 
over  the  hull  surface  to  obtain  the  thrust-deduction  force.  Since  the  interaction  force  is  to  be 
obtained  directly  from  relations  between  singularities,  the  main  problem  is  to  find  appropriate 
generating  singularities  for  hull  and  propellers.  (Essential  parts  of  the  problem  are  determina- 
tion of  the  perturbation  velocities  due  to  a hull  in  a uniform  potential  flow  and  determina- 
tion of  the  radial  distribution  of  propeller  thrust.  Specifically,  these  essential  parts  enter  into 
the  solution  of  the  interaction  or  thrust-deduction  force,  using  the  Lagally  theorem,  by  pro- 
viding the  required  hull-disturbance  velocity  and  propeller-sink  strengths. 

The  principal  assumptions  and  limitations  that  are  involved  in  the  analysis  and  treatment 
of  the  problem  are  summarized  as  follows: 

1.  Interference  effects  which  involve  changes  in  boundary  conditions  and  lead  to  itera- 
tions are  assumed  to  be  small4*13  and  are  not  considered. 

2.  Tlie  effects  of  various  singularities  can  be  combined,  lixact  linear  superposition  of 
flow's  is.  of  course,  limited  to  those  which  satisfy  the  Laplace  equation.  It  is  assumed  that 
the  effect  of  the  propellers  on  the  hull  is  mainly  p icniiul  in  origin,  whereas  the  effect  of 
the  hull  or.  the  propellers  is  essentially  viscous  in  ongl  The  effect  on  the  propellers  is 
treated,  however,  as  a potential  problem  through  the  use  or  the  circumferential  average  of 
the  total  wake  at  each  propeller  radius. 

3.  It  is  assumed  that  the  actual  propeller  blades  can  be  vplaced  by  lifting  lines  where 
the  induced  velocity  is  based  on  the  actual  radial  thrust  distributor,.  For  more  than  three 
blades,  calculations  show  the  circumferential  variation  (fluctuation)  ol  ‘minced  velocity  from 
lifting-line  theory  Iras  a negligible  effect  on  tire  thrust  deduction.  Lifting-lme  theory  is  used 
strictly  as  a matter  of  convenience  to  obtain  the  required  steady  sink  disk  strength  which  is 
radially  variable.  This  mathematical  propeller  model  gives  the  propeller  induced  velocity  as 
a function  of  propeller  thrust  in  a viscous  fluid  and  radius,  but  does  not  conside.  , opeller 
blade  thickness.  A check  of  the  effect  of  propeller  blade  thickness  was  marie  and  was  found 
to  be  small  for  the  propellers  considered. 

4.  It  is  assumed  for  this  investigation  that  the  separation  point  on  the  hull  is  not 
changed  by  the  propellers.  In  principle  an  iteration  procedure  and  boundary-layer  theory 
can  be  used  to  consider  changes  in  the  separation  point. 

5.  Tlie  theory  presented  is  applicable,  strictly  speaking,  only  to  ships  without  significant 
wavemaking  unless  the  potential  part  considers  the  free  surface.  To  make  the  theory  more 
general,  an  influence  coefficient  is  used  to  correct  the  potential  wake  for  the  ship  wave  system. 
In  the  present  report,  the  influence  coefficient,  also,  contained  some  geometrical  effects. 
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(>.  At  off-design  propeller  loadings  the  propeller  singularity  strengths  are  assumed  to  be 
proportional  to  (I  + 1/2  C’TI|).  This  assumption  restricts  the  analysis  to  eases  for  which  the 
relative  loading  difference  between  the  two  propellers  is  not  great. 


SOLUTION 

THRUST  DEDUCTION  AS  A LAGALLY  FORCE 

The  Lagally  steady-motion  equation16 


oo0 

'q  s^rr 

4rrr 


(1) 


connects  the  force  F(-  on  a source  at  some  point  in  the  How  field  with  the  output  Q and 
the  (low  velocity  Q(l/4rrr  at  the  location  of  the  source  Q.  where  r is  the  distance  from  a 
source  Q to  a source  Q.  For  contrarotating  propellers  Equation  (1)  can  be  rewritten  in  the 
more  general  form 


I'o  = PQ|  I V(w  )j  + (Uaj),  | for  forward  propeller 


(2) 


F0  = - p02  IV(w((),  + (Uaj)2 1 for  aft  propeller 

where  w is  the  hull  nondimensional-disturhing  velocity  and  U. , is  the  propeller  induced 
interference  velocity  and  at  each  conjugate  propeller  can  vary  spatially  but  not  temporally. 
The  hull-disturbing  velocity  wp  is  the  potential  flow  |>n«urbation  velocity  caused  by  the  hull 
moving  at  velocity  V.  The  hull  can  be  considered  as  generated  by  a surface  distribution  of 
sources  and  the  interaction  force  of  a contrarotating  propeller-hull  system  is  obtained,  by 
reciprocity,  from  the  hull-disturbing  velocity 


.-if 

P V i 


qm  (S)  ds 


V J A 1 

v $ 4»r 

at  the  contrarotating  propellers.  Each  propeller  is  represented  by  a sink  disk  of  strength 

<V  dA0 


l 


where  qm  = strength  of  a surface  source 
S = surface  area  of  the  hull  and 
A0  = propeller  disk  area 

Expressions  for  the  propeller  interference-axial  velocity  Uaj  will  be  given  later. 
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Equation  (2)  can  be  written  in  terms  of  wp  and  Uaj  at  each  propeller  disk  (of  the 
contrarotating  pair)  and  distributed  surface  sinks  of  strength  qm  = U . The  total  drag 

in  3 s 


augmentation  coefficient  CD  is  then  given  by 


2irpVI 2R,2 


f , <Uas>i 

= J [<Vi+  ~ — xdx 


/m!  t r,  w ,u^i  ,v‘A  . 

IV  J L Wp)!  “v  "J  ~ xdx 


■r 


(Vi 


(«..).  . j‘  r«U 


xdx  + 


xdx 


©' ) 


(Uas)2 

(Wp)2  -y~  Xdx 


(3) 


©'  f P^] 


xdx 


= Ij  + I,2  + I2  + I?1 

where  the  variables  within  the  integrals  arc  understood  io  be  average  values  in  the  circum- 
ferential direction. 

Interaction  force  has  been  presented  as  a drag-augmentation  coefficient  CD.  The  thrust- 
deduction  coefficient  is  related  to  the  drag-augmentation  coefficient  by 


t=tf  + t,  +t2=tf  + 


4<CDl+CD2> 


4Cr 


vThs 

where  tf  is  the  frictional  thrust  deduction.  For  tf  = 0 


"tf+  C 


Ths 


(V 


i 

t 


I 


and 

t 


'•'D  1 V‘D 

Some  comments  will  be  made  later  concerning  the  frictional  thrust  deduction  coefficient  tf. 
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The  interference  force  integrals  I,2  and  I21  in  Equation  (3)  cancel  since  only  the  force 
on  the  body  is  wanted.  The  simple  case  sketched  below  is  illustrative.  Assume  all  perturba- 
tions are  caused  by  three  point  sources  (representing:  body.  Propeller  1 , and  Propeller  2) 
enclosed  by  spheres  of  infinitesimal  radius  and  located  coaxially.  On  the  axis,  distance  r 
between  singularities  is  indicated  by  bar  (-).  From  the  system  as  sketched  we  can  write  the 
following  force  equations: 


VELOCITY  v INDUCED  AT  Q2  BY 
ALL  SINGULARITIES  EXCEPT  0, 


It  can  be  seen  from  liquation  (8)  that  the  interference  terms  Q,  Q2  on  the  .ight-hand  side 
cancel  since  the  distance  rjj  = r,^. 

INDUCED  VELOCITY  IN  CONTRAROTATING  PROPELLERS 

A design  procedure  for  contrarotating  propellers  which  is  based  on  Lcrbs’  theory  has 
been  presented  by  Morgan.17  Since  the  theory  and  calculation  of  wake-adapted  contraro- 
tating propellers  are  not  within  the  scope  of  this  report,18  the  procedure  used  here  for 
obtaining  the  propeller  source  strength  qm/V  = Uas/V  will  be  discussed  only  briefly.  To 
derive  the  thrust  deduction  from  the  Lagally  theorem,  it  is  important  to  emphasize  that  the 
required  singularity  strength  (far  downstream)  qm/V  is  equal  to  twice*  the  circumferential 
average  of  the  axial  component  of  the  propeller  self-induced  velocity  ;*t  the  propeller  disk.19 
It  was  shown  in  the  previous  section  that  the  propeller-induced  interference  terms  (U3i),  /V 
and  (Uaj)2/V  in  Equation  (3)  do  not  contribute  to  the  thrust  deduction. 

l7Morgan,  W.B.,  "The  Design  of  Counterrotating  Propellers  Using  Lcrbs’  Theory,"  Transactions  SNAMK,  Vol.  68  (I960). 

18 Lcrbs,  H.W.,  "Contra-Rotating  Optimum  Propellers  Operating  in  a Radially  Non-Uniform  Wake."  David  Taylor  Model 
Basin  Report  941  (May  1955). 

!9Wcinblum,  G.,  "The  Thrust  Deduction."  American  Society  of  Naval  Engineers,  Vol.  63  (1951). 

•The  actual  numerical  values  of  the  self-induced  axial  velocity  at  the  lifting  line  of  each  propeller  were  multiplied  by  two 
to  obtain  the  final  (far  downstream)  velocity  U^/V. 
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Lcrbs’  theory  for  either  lightly  or  moderately  loaded  contrarotating  propellers  considers 
propellers  having  finite  hub,  finite  number  of  blades,  and  an  arbitrary  radial  load  distribution. 
His  theory  is  based  on  potential  flow  in  which  viscous  effects  are  neglected;  however,  the 
drag  forces  on  the  blades  which  exist  in  a real  fluid  are  introduced  later  as  a correction  to 
the  potential  theory.  The  inflow  velocities  required  for  the  calculation  of  a wake-adapted 
propeller  are  based  on  the  total  wake.  The  nominal  inflow  Vx/V  is  corrected  by  V#/Vy  , 
which  is  the  ratio  of  the  effective  velocity  to  the  volume  mean  velocity.  The  volume  mean 
velocity  ratio  is  defined  as 


1 


xh 


xdx 


where  V /V  is  the  circumferential  average  of  the  local  velocity  at  each  radius. 

Calculations  for  moderately  loaded  contrarotating  (CR)  propellers  using  Lerbs’  induction- 
factor  method  are  programmed  for  a high-speed  digital  computer  at  this  Center.*  A 
numerical  solution  which  includes  the  total  propeller-induced  velocity  at  the  lifting  lines 
(see  Figure  1 for  velocity  component  diagram)  for  the  propeller  operating  in  a real  flow 
is  obtained  from  the  program.  Because  the  propeller  self-induced  velocity  Uas  is  needed  in  the 
calculation  of  the  thrust  deduction  and  only  the  propeller  total  induced  velocity  UaT  is 
included  in  the  computer  output  (Uas  not  convenient  to  retrieve)  Uas  must  be  calculated 
from  the  computer  output  by: 


(Ua4),  = (UaT)j  - (Uaj),  for  forward  propeller 


(Uas),  =fUa.,)2  (Uaj),  for  aft  propeller 

where  from  Morgan17 

<u3i>.  =^>2  <‘;>2  n (r,)2i 


(9) 


c;< 


J 2\  (U,/V) 
te,>l  *^>2 


= factor  for  obtaining  the  circumferential  average  of  interference  velocities 

= distance  factor  for  obtaining  the  effect  of  axial  distance  on  interference 
velocities 


* Reprogramming  of  die  present  OR  propeller  unnpntcr  program  is  currently  underway.  The  so-called  "equivalent 
propeller  com  opt”  ssill  no  longer  lie  needed  and  the  self-induced  velocity  |J  will  be  provided  in  the  output.  Consequently, 
ills-  method  used  to  obtain  ll  in  the  present  investigation  will  not  lie  necessary  in  the  future. 


8 


=r  ^ * gse^^T? 


*,W 


II 


G*  = total  nondimensional  circulation  at  each  radius  of  one  propeller  (ZP/ffDV) 

x = nondimensional  radius 

U(  = tangential  velocity 

D = propeller  diameter 

V = ship  speed 

Z = number  of  blades,  and 

P = circulation  of  each  section 


Tire  computer  program  uses  an  iterative-type  method  of  solution  to  determine  the  induced 
velocities.  Before  proceeding  with  the  design  of  two  actual  propellers  with  a specified  axial 
spacing  an  “equivalent"  propeller17  which  produces  one  half  the  total  thrust  is  introduced. 
The  expressions  for  the  self-induced  axial  velocities,  appearing  in  liquation  (9).  of  each 
actual  CR  propeller  in  terms  of  the  induced  velocity  U.  from  the  equivalent  propeller  are: 


(UaS>,  = <Ua>. 

(Uas),  =(!  +5  )(1  +F)(Ua)2 


GO) 


where  6 is  the  average  contraction  ratio  of  the  slipstream  at  the  aft  propeller  and  f is  the 
average  circulation  factor  at  the  aft  propeller.  According  to  Morgan17  both  6 and  f are 
small  in  absolute  value.  In  the  present  investigation  they  are  assumed  to  be  zero. 

Lerbs’  distance  factor18  ga  has  been  calculated,  for  the  present  investigation,  to  cover 
a higher  range  of  propeller  spacing  d/R  and  cunes  of  g versus  d/R  are  shown  by  Figure  2. 

M 

Recent  design  experience  at  the  Center  has  indicated  a preference  for  these  factors  over 
those  presented  by  Morgan.17 


THRUST  DEDUCTION  AT  OFF-DESIGN  LOADING 

The  minimum  horsepower  for  a ship  propelled  by  CR  propellers  is  a function  of  the 
propeller  efficiency  and  the  hull  efficiency.  Equal  thrust  or  equal  torque  on  each  propeller 
does  not  necessarily  give  the  minimum  shaft  horsepower  since  the  minimum  thrust  deduction 
docs  not  usually  occur  when  either  the  thrusts  or  torques  are  equal.  Therefore,  to  obtain 
minimum  horsepower,  it  is  necessary  to  investigate  the  effect  of  unequal  thrust  loadings  of 
the  propellers  on  the  thrust  deduction.  To  calculate  the  thrust  deductions  at  off-design 
loadings,  it  is  assumed  that  the  propeller  self-induced  axial  velocity  Uas  is  proportional  to 
the  propeller  thrust  coefficient  CTh. 

For  a particular  propeller  the  average  factor  fa  in  liquation  (9)  is  essentially  constant 
with  a change  in  thrust  only.  Thus,  for  a specified  CR  propeller  spacing,  i.e..  distance 
factor  ga.  we  observe  that  (Uai), = % constant  and  (U3 s )2 /f Ua% ),  = k,  **  constant 

and  the  off-design  relations  are: 
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2T,'  (Uas),  + k,  (Uas)2  2T2' = (Total  Axial), 

2T,'  <Uas>2  + MUas>I  2T1 ' = <Total  Axia»2 
At  design  condition  (assumed) 


The  necessity  for  utilizing  an  off-design  approximation  lies  in  the  fact  that  the  present  CR 
propeller  design  computer  program  (lifting-line  computation)  considers  only  equal  torque 
balance,  i.e..  T, ' - f,'. 

With  propeller  off-design  singularity  strength  taken  as  proportional  to  T'  in  Equation  (3), 
Equation  (4)  would  become 

8oyc  +t2'c) 

, = tf  + — (4a) 

’“Ths 

and  for  = 0 


where  the  drag  augmentation  coefficients  C,,  and  C,,^  are  at  design  condition. 

THRUST  DEDUCTION  AT  OFF  DESIGN  SPACING 

The  spacing  between  propellers  and  between  propeller  and  hull  can  also  affect  the 
minimum  shaft  horsepower.  As  soon  as  the  position  of  either  CR  propeller  is  changed,  all 
perturbation  velocities  must  be  recalculated  at  both  propellers.  For  off-design  spacing,  in 
this  investigation,  the  forward  propeller  location  remained  fixed  and  the  aft  propeller  was 
moved  downstream.  Large  propeller  spacing  (defined  by  axial  distance  between  propeller 
planes  as  a fraction  of  forward  propeller  radius.  d/R, ) requires  the  physical  removal  of  the 
rudder  for  the  typical  single-screw  ship.  Therefore,  thrust  deduction  was  determined 
mathematically  with  the  rudder  removed  when  the  parameter  d/R,  was  varied. 

To  approximate  the  change  in  the  hull-potential  wake  fraction  wp  with  distance  down- 
stream. a method  based  on  the  variation  of  thrust  deduction  with  propeller  axial  position 
behind  a body  of  revolution,  as  derived  by  the  author.4  was  used  together  with  the  exact 
functional  relation  between  the  potential  wake  fraction  ami  thrust-deduction  coefficient  for 
uniform  flow.  Tlie  relation  can  be  expressed  in  the  following  convenient  form: 
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Consider  the  fact  that  for  the  same  propeller  average  total  thrust,  we  have  the  same 
average  total  induced  velocity.  By  a judicious  choice  of  d/Rj  some  reasonable  assumptions 
concerning  off-design  performance  can  be  made  as  follows:  (1)  If  we  do  not  make  too  large 
a change  in  spacing,  the  assumption  is  reasonable  that  the  radial  distribution  of  thrust  for 
each  CR  propeller  remains  approximately  the  same.  With  this  assumption  we  calculate  new 
self-induced  velocities  Uas  at  the  new  location.  (2)  If  the  aft  propeller  is  positioned  down- 
stream where  wp  ^ 0 (about  d/R,  = 10.8  for  the  MAR  AD  Containers!#)  then  t,  = 0 and 
t ^ t, . It  is  known  that  T,  will  change  at  the  same  rpm  but  the  ratio  D/T  remains  almost 
constant  for  nonconfigurational  changes  in  propeller  loading. 


THRUST  DEDUCTION  DUE  TO  RUDDER 

A simple  relation  between  wp  and  tp  was  given  (liquation  (ID)  with  the  assumption 
of  uniform  flow.  Several  simple  and  reasonable  approximations  to  the  disturbing  velocity 


w induced  by  a rudder  may  be  made.  For  example,  a point  source  in  a two-dimensional 
uniform  flow  (represents  a rudder  of  in  Unite  span  and  chord)  may  be  an  adequate  model 
for  the  present  investigation  inasmuch  as  only  the  flow  field  ahead  of  the  rudder  is  needed 
and  points  beyond  the  propeller  slipstream  lack  a point  of  application.  Isowake  lines  can 
be  constructed  at  the  propeller  plane  and  these  integrated  circumferentially  at  several  radii 
to  provide  an  average  radial  distribution  of  wp . 

Harvald20  has  constructed  a two-dimensional  flow  model*  by  a conformal  transforma- 
tion technique  that  produces  a more  typical  rudder  section  shape  than  a two-dimensional 
half  body.  His  Diagram  89,  Chapter  V.  was  used  to  evaluate  wp  for  the  MARA1)  Container- 
ship  rudder,  liquation  (II)  was  then  used  to  calculate  the  thrust  deduction  due  to  the 
nulder. 

COMPUTATIONAL  EXAMPLE 

MARAD  CONTAINERSHIP 

Figure;.  3 and  4 show,  respectively,  the  lines  and  stem  view  of  the  Center  Model  5218, 
representing  a 25.5-knot  containeiship  equipped  with  a CR  propeller  stem  arrangement 

■^Harvatd,  S.A..  “Wake  of  Merchant  Ships.”  Hie  Ifanish  Technical  Press,  Copenhagen  (1950). 

Mlarvxld  also  examined  the  effect  of  a convergent  flow  field  and  it  appears  that  in  general  an  assumption  of  a free 
rudder  in  a parallel  How  would  be  adequate  for  the  purpose  of  the  present  investigation. 
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and  a single  rudder.  Tins  con  tain. 'rsilip  has  been  investigated  extensively  by  Strom-Tejsen14 
with  model  tests  of  several  desi  ^n  variations  including  twin  rudder?.,  single  screw,  and  over- 
lapping propellers.  Ship  and  model  data  for  this  model  are  given  in  Table  1 . A drawing  and 
pertinent  geometrical  data  appear  in  Figure  5 for  the  design  CR  propellers,  Numbers  4458 
and  4459.  Surface-sink  strengths  of  these  propellers  as  a function  of  propeller  radius  for  the 
design  condition  (equal  thrust)  and  at  an  off-design  spacing  d/R,  = 1.5  are  given  in  Figures 
6 and  7,  respectively.  Table  2 gives  predicted  propulsion  results  at  25.5  knots  for  model 
experiments  for  the  ship  equipped  with  contrarotating  propellers.  The  experiments  conducted 
with  Model  5218  included  propulsion  data  for  the  conditions:  with,  and  without  rudder, 
off-design  thrust  loading  between  CR  propellers,  and  comparative  single-screw  results.  Titus, 
it  seemed  to  be  an  excellent  vehicle  for  performing  a theoretical  analysis  of  the  thrust 
deduction  for  a high-speed  single  screw  ship. 

INFLUENCE  COEFFICIENT  FOR  POTENTIAL  WAKE 


To  effect  a major  saving  in  time  and  money  the  potential-wake  fraction  wp  was  obtained 
for  the  subject  containershsp  by  application  of  an  influence  coefficient.  The  influence 
coefficient  would  not  generally  be  needed.  In  essence  this  coefficient  was  determined 
numerically  from  propulsion  data14  (Table  2 of  this  report)  and  the  experimental  radial 
distribution  of  wp  reported  for  the  single-screw  surface  ship  SIMON  BOLIVAR.10 

Introducing  the  influence  coefficient  e into  F.quation  (4).  we  can  write 

4cC, 


^experimental  V 4" 


I) 


■ 1 hs 


(4b) 


where  eC„  = efl , + I2)  by  Equation  (3)  and  the  integrals  I,  and  I2  have  been  computed  for 
the  MARAI)  Containership  using  a wp  distribution  from  SIMON  BOLIVAR  already  adjusted 
Equation  (1 1)  for  differences  in  propeller  diameter  and  spacing  relative  to  the  hull.  Let 

tf  = 0.015  from  Beveridge10 


(t)CJ 


M’ 


= 0.188  from  Table  2 


C,hs  = 0 6924  from  Strom-Tejsen14  (actual  test  value  for  CR  set)  and 
= 0.0240  from  Equation  (3) 

Equation  (4b)  solved  for  c is 

(1hx  lft,e\p  ~ M 
C ' 4C„ 


= 1.25 


(12) 


The  influence  coefficient  c maip’.y  corrects  the  assumed  wp  distribution  for  differences  in 
the  free-wave  system  and  geometry  between  SIMON  BOLIVAR  and  the  MARAD  Container- 
ship.  Tlie  effect  of  geometry  appears  to  be  dominant.  Although  other  unaccountable  effects 
may  be  lumped  into  e it  is  believed  to  operate  chiefly  on  wp  (as  intended)  by  a comparison 
with  wake  data  of  Harvald.21  and  Nowacki  and  Sharma.22  Final  curves  of  w(),  and  wp, 
versus  x for  the  containership  are  presented  in  Figures  8 and  9 at  the  design  CR  propeller 
spacing  d/R(  = 0.52  and  at  off-design  spacing  d/R,  = 1.50. 

A check  of  the  Harvald  data  reveals  w = 0.22  « wp  for  a small  wave  wake  while  the 
present  containership  potential  wake  (Figure  8)  shows  wp  = 0.21  at  0.7  propeller  radius. 
Another  check  on  wp  for  the  containership  can  be  obtained  from  the  relation  of  t and  wp 
for  uniform  flow  (Equation  (1 1 )).  With  C1})  = 1.214  (computed  for  1 - w()  = 0.76)  and 
t = 0.188  - 0.015  = 0.173.  we  have 

W„  - f (I  ♦v'l  +CTh)  =0.22 

A small  wave-wake  component  w^  is  indicated  at  design  Froude  number.  F^  = 0.272.  for 
the  containership  and  some  substance  is  given  to  this  by  the  results  of  Nowacki  and  Sharma, 
which  are  given  in  Figua*  10.  for  a somewhat  idealized  mathematical  hull  having  = 0.64 
and  Cx  = C'p  = Cw  = 0.8. 

NUMERICAL  RESULTS 

Although  the  frictional  component  of  the  thrust-deduction  coefficient  tf  was  introduced 
earlier,  it  was  assumed  equal  to  zero  in  all  numerical  computations  involving  ratios  of  thrust 
deduction  (e.g..  t,  /t  and  t,/t).  This  assumption  avoids  any  superficial  division  such  as  (tf), 
and  (tf),  for  an  already  small  quantity.11,13  However,  ah  absolute  values  of  t are  based  on 
t,  = 0.01 5. 13  A most  important  observation  is  made  at  this  point  concerning  the  influence 
coefficient  c that  was  used  in  determining  potential  wake:  to  wit.  c cancels  if  we  deal  with 
thrust  deduction  ratios  in  our  analysis.  In  a very  real  sense  percentage  or  fractional  changes 
in  thrust  deduction  with  parametric  variations  are  what  we  want  anyway. 

RUDDER  CONTRIBUTION  AND  SPACING  BETWEEN  PROPELLERS 

As  described  previously  the  Harvald2”  diagram  showing  the  dependence  of  wp  on  the 
dimensions  and  position  of  the  rudder  (in  parallel  flow)  may  be  used  in  conjunction  with 


‘ HarvalJ,  S,.'.,  "Wake  anil  Tbru\t  IWuction  at  l-.Mrcmr  Propeller  l.oadinpv"  Sweiliili  Stale  Shipbnililin^  l:\perimental 
tank  Publication  \<»  61.  Ooteborp  (1967). 

“Nowacki.  II.  anil  S.l).  Sharm.i.  ”1 1-  c-Stirfuce  Iffccts  in  Hull  Propeller  Interaetion.”  The  Universil>  of  Miehip.in 
(’ollcj’c  of  l.npincerinp  Report  112  (Sep  1971). 
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liquation  (11)  to  estimate  the  thrust  deduction  due  only  to  the  rudder  For  the  MARAD 
Ccntainership  the  following  rudder  input  has  been  used: 

t*/C  = 0.43 
D/«  = 2.6 
a/U  = 0.45 

where  t*  = maximum  section  thickness  of  the  rudder  (symbol  t*  not  to  be  confused  with 
thrust  deduction  coefficient), 

D = propeller  diameter 

a = axial  clearance  between  the  rudder  leading  edge  and  a representative  propeller 
plane.  Note:  Propeller  plane  defined  at  0.71)  when  rake  is  present  and 
(’  = rudder  chord  length 

For  the  subject  rudder  and  C R propeller  set  the  midplane  was  chosen  in  computing  (a),  and 
(G)  was  taken  as  the  rudder  chord  at  0.7Rj . A w()  of  0.05  results  and  (tp)  rudder  = 0.04 
was  computed  from  liquation  (II).  The  experimental  result  for  the  ccntainership  from  tests 
with  ami  without  rudder  gave 

At  = 0.056  Design  CR  propellers  (Table  2) 

At  = 0.029  Single  Screw.  Strom-Tejsen14 

Very  good  agreement  between  the  experimental  and  theoretical  results  for  the  thrust  deduc- 
tion due  .solely  to  the  rudder  is  revealed  when  it  is  realized  that  choosing  a midplane  location 
for  the  CR  propellers  is  tantamount  to  placing  the  propeller  disk  loading  closer  to  the  rudder 
compared  to  the  single  screw  location  and  moving  the  aft  CR  propeller  farther  ahead  of  the 
rudder.  Titus,  the  computed  incremental  value.  At  = (tp)ruJdl,r  = 0.04,  lies  between  the  two 
experimental  values. 

A simple  monotonic  variation  of  the  thrust  deduction  coefficient  with  axial  distance 
between  CR  propellers  was  found  as  expected  (see  Figure  1 1).  The  rudder  has  been  removed 
(physically  necessary  for  large  spacing)  by  subtracting  its  known  contribution  from  the  total 
thrust  deduction  since  it  had  been  determined  experimentally.  Selection  of  particular  discrete 
d/Rj  ratios  for  performing  thrust  deduction  computations  was  discussed  previously.  There 
is  nothing  striking  or  unusual  about  the  results  obtained  and  the  variation  for  the  MARAD  CR 
propellers  is  similar  to  the  curve  (shown  in  Figure  1 1 for  comparison)  derived  for  a bare-hull, 
single-screw  submarine  form.4 

Although  the  present  work  was  limited  to  an  investigation  of  thrust  deduction,  some 
more  far-reaclung  comments  seem  desirable  at  this  point.  It  is  generally  accepted  that  axial 
spacing  has  little  effect  on  the  open-water  efficiency  of  CR  propellers,  provided  they  are 
operated  at  their  design  spacing.  Limited  space  available  usually  dictates  that  propeller 
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spacing  be  kept  to  a minimum.  Meeker23  has  investigated  CR  propeller  performance  in  uni- 
form flow  for  distance  ratios  from  d/Rj*  0.2  to  0.8.  Purely  from  a hydrodynamic  standpoint, 
the  concept  that  behind  a hull  the  potential  flow  field  decays  at  a greater  rate  than  does  the 
frictional  flow  field,  is  a basis  for  speculating  that  an  optimum  propeller  spacing  might  exist 
in  the  wake-adapted  case.  As  an  example  it  has  been  shown  for  a submerged  body  of  revo- 
lution4 that  an  optimum  (based  on  maximum  propulsive  coefficient)  downstream  location 
does  occur  for  a single  screw.  Based  on  these  statements  a spacing  for  optimizing  the  pro- 
pulsive coefficient  probably  exists  for  a CR  propeller  pair  installed  on  a submersible  and 
certain  surface  ships  where  free  surface  and  wave  effects  would  not  be  significant. 

Practical  application  of  these  ideas  may  be  feasible  because  submcrsibles  arc  presently 
configured  with  the  stem  control  surfaces  forward  of  the  propeller,  and  a twin-rudder  arrange- 
ment for  a surface  ship  would  permit  the  aft  propeller  of  a CR  set  to  be  positioned  farther 
downstream.  In  retrospect,  investigative  effort  along  these  lines  seems  warranted. 

THRUST  BALANCE 

Table  3 gives  calculated  thrust  deduction  results  at  several  discrete  values  of  relative 
thrust  loading  between  the  forward  and  aft  propellers  designed  for  the  M ARAD  Containcrship. 
For  equal  thrust  loading  T,  /T2  = 1 .0  the  division  of  the  thrust  deduction  is  seen  to  be  73 
percent  of  the  total  produced  by  the  forward  propeller  and  27  percent  produced  by  the  aft 
propeller.  Tins  result  may  be  important  in  the  propeller  design  problem  since  an  assumed 
50-50  distribution  of  the  thrust  deduction  would  thus  produce  approximately  a 5 percent 
difference  in  the  division  of  useful  propeller  thrust  between  propellers. 

As  mentioned  earlier  a design  condition  T,/T2  = 1.0  was  assumed  in  the  theoretical 
analysis.  However,  the  experimental  results  for  the  MARAD  Containcrship  reported  by 
Strom-Tejsen14  revealed  that  at  equal  rpm  a value  T,  /T2  = 0.95  was  obtained  in  contrast  to 
equal  thrusts  for  which  they  were  designed.  In  Figure  12  the  computed  (1  - t)  values  of 
Table  3 have  been  plotted  as  fractions  of  the  design  value  with  an  incremental  shift  of  0.05 
on  T,/T2  to  give  a unity  thrust  deduction  ratio  at  T,/T2  * 0.95.  In  addition  two  other 
experimental  points  are  shown  corresponding  to  the  rpm  ratios  of  0.95  and  1.05  in  Table  2. 
The  (1  - t)  ratio  curve  of  Figure  12  shows  about  an  8 percent  decrease  when  the  propeller 
thrust  ratio  T,/T2  is  increased  from  about  0.4  to  2.0.  The  two  off-design  experimental 
points  do  not  lie  on  the  theoretical  curve.  This  is  believed  to  be  due  to  either  experimental 
test  accuracy,  or  the  rudder  effect,  or  to  both.  In  percent  the  differences  are  small,  but  it 
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decker,  R.  and  N.A.  McDonald,  "The  Axial  Spacing  and  Optimum  Diameter  of  Counterrotating  Propeller*,”  David 
Taylor  Model  Baxin  Report  1 342  (Feb  I960). 
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should  be  pointed  out  that  the  rudder  effect  on  thrust  deduction  was  treated  theoretically 
as  a constant  disturbance  factor  in  tin.  singularity  system,  while  in  physical  reality  one  might 
expect  a thrust-loading  condition  that  would  produce  a minimum  thrust  deduction  with  a 
rudder.  This  problem  requires  further  investigation.  For  the  present  case  experimental  data 
at  moa*  extreme  propeller  thrust  loadings  would  have  been  informative.  The  earlier  discussioi 
concerning  the  rudder  contribution  to  the  thrust  deduction  and  its  overall  effect  on  CR 
propulsion  (Table  2)  seem  to  show  that  the  single  centerline  rudder  offers  no  advantage  in 
this  propulsion  system. 

Also,  depicted  in  Figure  12  is  the  theoretical  variation  of  the  separate  thrust-deduction 
coefficients  versus  T,/T,  for  the  forward  and  aft  propellers.  Of  special  interest  is  the  inter- 
section of  the  curves  for  t(  = t,  which  shows  that  a thrust-loading  ratio  of  about  0.4  is  needed 
to  accomplish  equal  division  of  the  thrust  deduction  in  this  CR  propeller  set.  The  curves 
have  been  extended  to  the  limiting  values  as  T,  -*  0.  Reverse  limiting  values  are  obtained  at 
the  upper  limit  T,  -*  A calculation  shows  that  t,/t  = 0.98  and  t,/t  * 0.02  at  T,/T2  = 20. 
These  statements  concerning  extreme  values  of  T|  / T,  have  only  academic  meaning,  since 
they  are  not  practical.  However,  the  required  linearity  (Assumption  (>.  Statement  of  Problem) 
is  preserved  because  the  total-thrust  coefficient 

r,  + t, 

C„1S  , — = 0.6924 

T*Vl  V' 

is  constant. 

Determination  of  the  separate  thrust  deduction  coefficients  t,  and  t,  for  the  CR 
propellers  was  essential  to  finding  the  total  thrust  deduction  of  the  CR  propeller  system  and 
this  point  is  emphasi/.ed.  The  variation  of  t(  and  t,  with  the  thrust  ratio  T,/T,  may  be 
important  to  the  total  propulsion  problem.  As  discussed  in  the  section  on  THRUST 
DEDUCTION  AT  OFF-DHSlCiN  LOADING,  equal  thrust  or  equal  torque  may  not  necessarily 
give  minimum  shaft  horsepower.  For  iIk  containcrship.  Figure  12  shows  that  the  total  thrust 
deduction  is  reduced  b\  increasing  the  thrust  on  the  aft  propeller  relative  to  the  forward 
propeller.  It  may  be  possible  to  design  a C'R  propeller  pair  (without  centerline  rudder)  for 
unequal  thrust  producing  a lower  thrust  deduction  and  no  less  in  propeller  efficiency  com- 
pared to  an  equal  thrust  condition.  That  this  is  possible  is  not  reflected  by  the  horsepower 
given  in  Table  2 for  the  containcrship  when  the  rpm  ratio  n,/n,  was  varied.  However,  it  is 
noted  in  this  case  that  the  thrust  unbalance  is  relatively  small,  the  rudder  effect  is  unclear 
and  the  off-design  propeller  efficiency  may  not  be  the  best  efficiency  obtainable. 
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CONCLUDING  REMARKS 


A theoretical  method  was  developed  to  determine  the  steady  interaction  force,  thrust 
deduction,  of  contrarotating  propellers.  Numerical  values  of  this  force  are  needed  to  de- 
termine the  propeller  thrust  required  to  propel  a ship  at  design  conditions,  and  to  make 
parametric  investigations  of  propulsive  performance.  In  the  method  the  force  was  obtained 
as  a Lagally-type  force  which  has  its  basis  in  a principle  of  reciprocity  for  finding  forces 
between  bodies  in  a flow  Held.  Certain  interference  forces  were  shown  to  cancel.  The 
methodology  involved  a combination  of  potential  and  viscous  flow  effects. 

Application  of  the  method  was  made  to  secure  numerical  results  for  a set  of  CR  pro- 
pellers designed  for  a MARAD  high-speed  containership.  Some  salient  features  and 
implications  concerning  the  thrust  deduction  of  the  containership  were: 

1 . Determination  of  the  separate  thrust  deduction  for  the  forward  and  aft  propeller 
ot  the  CR  pair  at  design  and  off-design  loading.  At  design  loading  73  percent  of  the  total 
thrust  deduction  was  produced  by  the  forward  propeller  and  27  percent  by  the  aft.  The 
total  thrust  deduction  was  in  exact  agreement  with  the  experimental  result  because  of  the 
use  of  an  influence  coefficient  in  determining  the  final  potential  wake.  Indications  are 
that  the  distribution  ot  the  thrust  deduction  between  propellers  might  Ik*  important  with 
regard  to  minimum  shaft  horsepower.  However,  additional  investigation  will  be  required  to 
clarify  this  point. 

2.  Computation  of  the  rudder  contribution  shows  agreement  with  experiment  within 
better  than  two  percent  for  one  minus  the  thrust  deduction. 

3.  As  expected  a simple  monotonic  reduction  in  thrust  deduction  was  observed  when 
the  forward  propeller  location  remained  fixed  and  the  aft  propeller  was  moved  downstream. 
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Figure  10  - Wake  Components  for  an  Idealized  Hull  Form  (From  Reference  22) 


Thrust  Deduction  as  a Function  of  Relative  Thrust  Loading 
between  CR  Propellers  for  MARA!)  Containership 


TABLE  I SHIP  AND  MODEL  DATA  FOR  MODI-L  5218 
CONTRAROTATING  CONTAIN ERSIIIP  DESIGN 


APPENDAGES:  Propeller  Bossing  (No  Rudder  or  Bilge  Keels) 


DIMENSIONS 


length  (LWL),  ft 
Length  (LWP),  ft 
Beam  (Bx),  ft 
Draft  (H).  ft 
Displ,  tons 
Wetted  Surf,  srjft 
Design  v,  knots 


394.8  Aft  of  FP 


LCBlbr  = 395.5  Aft  of  FP 


LWL  COEFFICIENTS 


Ship 

Model 

CB 

0.555 

CWF 
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cP 
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780.0 
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I 

Cx 
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le/l 

0.52 
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3.408 

cw 

0.697 

Lp/L 

0.00 

30.0 

0.989 

CPF 

0.57 

Lr/L 

0.48 

38520SW 

1.342FW 

CPA 

0.62 

L/Bx 

7.53 

83920 

91.7 

CPE 

0.58 

Bx/H 

3.45 

25.5 

4.63 

CPR 

0.60 

A/(OilF 

1 81.1 

u.ou  1545 

^ 0.74  f 009" 

fi  ne  * n nc. 
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X = 30.334  VVT^7l  = 0.914 
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